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ABSTRACT 


The Ships Parts Control Center (SPCC) Uniform Inve 
Control Program (UICP) wholesale replenishment model for lH 
cognizance symbol (consumable) material is an order quantity- 
reorder level or (Q,r) model. A stocked item's order quantity 
and reorder level are established in large part by the unit 
price and procurement lead time forecasted for it. When a 
replenishment is needed, the order quantity is specified and 
the procurement officer requests bids from vendors. These 
bids include both a unit price and production lead time. This 
thesis analyzes the influence of different bids with differ- 
ent unit price and different lead time on the future optimum 
total annual cost of stocking the item as computed by the 
UICP model. Based on this analysis, a simple technique to 
evaluate those bids is developed and steps to implement this 


technique are suggested. 
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۳ THE PROBLEM 

The current practice at the Navy's Inventory Control 
Points (ICPs) is to treat the wholesale level inventory manage- 
ment of an item and the procurement of replenishment stock as 
Separate functions even though they are, in fact, key elements 
of the same single supply system and share the goal of maximum 
fleet support within annual budget constraints. The Uniform inven- 
Ex sGontrol Point (UICP) inventory models are used to determine 
the reorder level and order quantity based on data from recent 
past procurements. When the inventory position of an item 
drops below the computed reorder point, an order is sent to 
the procurement ο 0 to buy the computed order quantity. 
The procurement department then solicits bids from potential 
vendors who are interested in filling the order. The vendor 
selected is usually the one who has bid the lowest unit price 
and can deliver within an ICP established desired delivery 
date. Some time later the UICP inventory models receive the 
new price and production lead time values of the vendor winning 
the contract. A new reorder point and order quantity are then 
computed and form the basis for the next procurement replenish- 
ment stock. 

If the information from the vendors' bids could be incor- 


porated in the inventory models at the time that the bids are 


being evaluated rather than after a bid has been accepted, a 
better evaluation of the bids should result. 126 Ότο. 
accepted should be the one which provides the lowest future 


average annual inventory costs. 


B. BACKGROUND 

For each item managed, the Navy's UICP models compute the 
values of the reorder point and order quantity which minimize 
the average annual variable costs of ordering, holding); و‎ 
time-weighted backorders. These models are based on the 
traditional steady-state continuous review lot size-reorder 
point models for stochastic demands which assume a constant 
price and procurement lead time for each item. Under the 
assumption of a constant price, the average annual total pro- 
curement cost is a constant value which is independent of the 
decision variables, order quantity and reorder point. As a 
conseguence, this cost term can be ignored. 

If, however, the procurement cost is known to be variable 
between vendors then the average annual procurement costs should 
be added to the average annual costs of managing the inventory 
of an item when a reprocurement is being considered since its 
magnitude is usually much larger than the average annual variable 
costs. The UICP models would then take on the form of stochas- 
tic price-break models. Solution techniques for solving such 
models are well known. Hadley and Whitin present these tech- 


niques in Chapters 2 and 4 of Reference 1. The UICP models were 


originally programmed to consider such price-breaks. In recent 
years, however, that option has not been used because of 
limited computer capacity. 

Another dimension of the bidding process is that vendors will 
typically submit different estimates of production lead time. 
Fortunately, inclusion of variable production and therefore 
procurement lead times in the UICP model does not require any 
additional average annual cost terms. The impact of lead 
time on the optimization is concentrated in the determination 
of the reorder point. However, the determination of the reorder 
point also includes an item's price because the reorder point 
influences not only the expected number of shortages but also 
the expected inventory holding costs. 

Hadley and Whitin [Ref. 1] have developed a price-break 
model for considering the influence of both unit price and 
lead time. In this model they assume that the unit price is 
both a function of the order quantity and the lead time. They 
assume, in particular, that the unit price offered by a vendor 
will increase as lead time decreases in a stepwise fashion; 
i.e., the same unit cost applies over a range of lead time 
values as well as order quantity values. An algorithm is also 
presented for solving this problem. Unfortunately, the model 
differs significantly from that of the UICP; there is no back- 
Order term in the holding costs and the shortage costs are 
based on the expected number of backorders rather than time- 


weighted backorders. In addition, while their assumed 


relationship between price and lead time might be true for 
negotiations with a single vendor, it is not obvious m EE 
would be valid for all potential vendors. 

Rather than attempt to identify if some such relationship 
exists for all vendors or to assume it even if it doesn't, 
it seems more efficient to conduct a study of the impact on the 
UICP model of varying combinations of unit price and procurement 
lead time and then to develop a simple methodologv for comparing 
two or more bids. The result would provide a procurement 
manager with a simple technique which would enable him to 
integrate the inventory management and procurement activities 
in order to minimize the total expenditure of Navy dollars re- 


quired to stock consumable items. 


۳ 5 

The objective is to develop a management technique for use 
by procurement personnel which will permit evaluation of vendor 
bids for consumable items on the basis of their impact on 
total average annual inventory costs. This technique should 
be both simple and quick to use with minimal requirement for 


computer or calculator equipment, 


D. APPROACH 

The approach will be to examine the impact of varying com- 
binations of unit price and procurement lead time values on total 
average annual costs predicted by the UICP consumable model to 
determine what savings are possible and how to achieve them. 


The effect of changing either the unit price or the 7 


10 


will be examined by taking partial derivatives of the total 
cost equation. In addition, the relationship between these 
variables and their effect on predicted inventory costs in the 
UICP consumables model will be examined through the use of a 
computer program which will first duplicate the UICP(Q,r) 
solution process, then incrementally change the variables of 


interest and compute the resulting inventory costs. 


E. SCOPE AND LIMITATIONS 

The stock procurement process will be examined at the point 
where the buy quantity has been determined and vendor bids have 
been solicited but no contract awarded. No required delivery 
date will be assumed. Since current procurement procedures 
do not retain bid data, the variations in lead time and unit 
price between bids are not known, and therefore the thesis 
model will allow for a large range of possible values of each. 

The UICP inventory model for SPCC managed lH cognizance 
consumable material will be utilized as the basis for the 
thesis model. For simplicity in programming and to keep the 
scope to manageable size, only items having sufficient average 
demand quantities such that their lead time demand quantities 
can be assumed to be normally distributed are considered. 
Slower moving items with Poisson or negative binomial distri- 
butions of lead time demand can be similarly analyzed with 
appropriate changes in the sections of the computer program 
which calculate the reorder level and the expected number of 


backorders. As a matter of fact, most of the replenishment 


1: 


procurement which is done by the Navy is for the normally 
distributed items even though they represent only 4$ of the 


items managed by SPCC and 13$ managed by ASO [Ref. 2]. 


F. PREVIEW 

Chapter II will present a brief overview of the current 
UICP consumables procurement model to establish the basis for 
analyzing it in both Chapters II and III. Model assumptions, 
constraints, and the total expected annual variable costs (TVC) 
equation and its optimization methodology will be discussed. 
Mathematical analysis of the TVC equation by taking partial 
derivatives and examining its behavior graphically will also 
be presented. 

Chapter III will describe the computerized incremental 
analysis of TVC when it is subjected to changes in the unit 
price and procurement lead time. The results will be then 
analyzed by using the isocost technique. The results will be 
illustrated graphically and limitations on the values of the 
bid variables will be suggested. 

Chapter IV will present a technique for evaluating bids 
based on the analyses of Chapter III. 

Chapter V will provide a summary of the chapters, present 
conclusions regarding the value of the bid evaluation tech- 
nique, and recommend steps for implementation. Recommendations 


for further research on the technique will also be provided. 
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II. THE CURRENT UICP MODEL 


A. GENERAL DESCRIPTION OF THE MODEL 

The Navy's Uniform Inventory Control Program (UICP) 
wholesale consumables model, used to set inventory levels for 
SPCC managed 1H cog items, forms the basis for the model 
developed in this thesis. The model seeks "to minimize the 
total of variable order and holding costs subject to a con- 
straint on time weighted, essentiality-weighted requisitions 
short" [Ref. 3]. The average annual total variable cost (TVC) 
equation used contains three main terms: an ordering cost 
term, or average number of orders per year times the adminis- 
trative cost to place an order; a holding cost term, or the 
average number of units on hand at any random point in time 
multiplied by the cost to hold a unit in stock for a year; and 
a shortage cost term, consisting of the average number of 
requisitions backordered at any random point in time multiplied 
by the cost incurred by not filling a requisition for a year 
times the military essentiality or worth of the item. The 
average annual cost of the items procured (unit price multi- 
plied by average annual demand) is considered a fixed cost 
independent of the decision variables and is not considered 


in the model. 
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The following assumptions apply to the UICP model. 


These assumptions will also be used in developing the technique 


will be presented in Chapter IV. 


Steady state environment--The mean and standard 
deviation of the random variables, quarterly demand 
and procurement lead time, are assumed constant over 
all future time. 


No quantity price discount--The unit price is the same 
regardless of the number of units in an order. A 
price-break subroutine is contained in the UICP 
implementation but it is not used at present. 


Instantaneous reorder--Replenishment orders are placed 
immediately after the inventory position drops below 
the reorder level. Although this assumption is a 
practical impossibility, the actual time delay is com- 
pensated for by including the associated administrative 
lead time as part of the procurement 1 me 


The cost to hold one unit of stock is proportional WEN 
the unit price of the item (currently set at 23% of 
the unit price per year). 


The time-weighted cost of a backorder for an item can 
be accurately quantified for determining stockout 
costs. Although this value (lambda) is actually 
determined from budget and supply material availability 
(SMA) constraints, for computational and analysis 
purposes, lambda will be assumed to accurately 
represent actual stockout costs. 


The military worth (essentiality) of an item can be 
accurately quantified, as required for the determina- 
tion of stockout costs.  Essentiality is currently 
fixed at 0.5 for all items by SPCC. 


No interaction exists between items. Each item's 
order quantity and reorder point can be determined 
independently of other items. Similarly, total inven- 
tory costs for a group of items can be determined by 
adding the independently computed costs for each item. 


which 


(1) 


(2) 


(3) 


(E) 


(5) 


(7) 


C. TOTAL VARIABLE COST EQUATION 


The UICP total average annual variable cost equation is 


presented below, with the first term representing the order 


l4 


cost, the middle term the holding cost, and the last term 


the backorder cost. 


AJA + I-C[R + Ë - L-D + B(O,R)] + AE[B(Q, R)] 


total variable costs of one stocked item 
per year. 


expected or average number of units demanded 
per quarter; forecasted from historic demand 
quantities and trends. 


order quantity. 


administrative cost of a procurement action; 
equal to $380 for purchases under $8,000, 
$1,050 for negotiated contracts (over 
$8,000), and $1,080 for advertised contracts 
(over $8,000). 


reorder level (based on inventory position, 
NOE JUSE SCOCK ON hand). 


procurement lead time (mean value forecasted 
from past procurement actions). 


expected number of units backordered at any 
Fandom point in time (a function of Q and R). 


annual inventory holding cost rate, composed of 
storage, obsolescence, and opportunity costs 
as percentages of unit cost (equal to 0.23 

for consumable items). 


ποτε cost of the item. 


expected number of units demanded per customer 
requisition. 


Shortage cost of one requisition backordered 
for one year. Currently set at $1,500 for 
category A (formerly lHOl and lHO2 cog) 
items, $1,000 for category B (formerly 

lHO3 cog) items, and $500 for category C 
(formerly lHbb cog) items. 


15 


TVC = 


where: 


TVC 


B(Q,R) 


E = military essentiality of the r1tem;rcourrerr m 

8٣۳٣ 
2 = average number of procurement actions or 
inventory cycles per year. 
[R +5 -LD +B(Q,R)] = expected number of units in stock 
at any random point in time (average on-hand 
inventory level). 
B(Q,R) 


————- = expected number of requisitions on backorder 
at any random point in time. 


The formula for B(Q,R) is, from Hadley and Whitin [Ref. 4], 


00 


BOUM NE f (w -R) [F(x +Q;L) - F(x:L)]dx 


where 


F(x:L) = the distribution function of deman T: 
over lead time L. 


Hadley and Whitin reduce this formula to the general form of 
T 
B(Q,R) - g [8 (R) μου 
When demand over procurement lead time is assumed to be normal 


with a mean of DL and a standard deviation of o, 


SD S : 2 R-DL, _ 8 _۔م‎ R-DL 
B(R) = x lo + (R-DL) “] 4 ( 3 ) 3 (R DL) $ (=> ) 
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and 


2 RYO =DE 


0 


R+Q-DL 


7 ۰۰01 


8(R+0) = اچ‎ ) - J(R+Q-DL) 9 ( ) 
The notation » corresponds to the density function for the 
normal distribution and ð corresponds to its complementary 


Sumulative distribution function. 


D. OPTIMIZATION AND KEY VARIABLE RELATIONSHIPS 

As with other inventory models, the UICP cost equation is 
minimized by taking the partial derivatives of TVC with respect 
to the decision variables, Q and R, and setting them equal to 
zero. Unfortunately, the results are two complex equations in 
Q and R. The solutions derived by the Operations Researchgroup at 
the Navy Fleet Material Support C ۴ لت‎ ۹٦٥٦6 ۰۶۹۱٥۸٥569 


values, and parameters needed to solvethe formulas areprovided below: 


OF = ry CAD IG < the usual EOQ 
DIC 
— = * 
RISK = |] F(R*) DIC ΤΠΕ 


and, from the table of the cumulative density function for the 
Standard normal distribution, we can find the value of the 


normal deviate, Z. R is then computed from 


E ο τα, 


where o is the standard deviation of demand during procurement 


lead time. The formula for c depends on the Mark Code of the 


7 


item. Mark II and Mark IV items are both characterized by 
having quarterly demand of more than 5. The difference between 
these "marks" is a function of the product of the 5 m, 
demand and the unit cost; Mark IV items correspond to items 
for which this product is greater than 75, Mark II items have 
product values not exceeding 75. 

For Mark IV items the formula for o is 


o = 41.57 MAD? .L DOM MAD? p^ 


where 


ΜΑΡ. = Mean absolute deviation of quarterly 
demand for item; forecasted from historical 
demand data. 

MAD, < Mean absolute deviation of procurement lead 
time; forecasted from prior procurement 
actions. 


For Mark II items the formula is 


c = V1.57 MAD2.r + D2.L 


D 


The optimal values of Q and R are then subjected to the 


following constraints: 


1:5; 
O = - 


| 
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O99 for category C items; 
0.01 < RISK < 0.40 for category B items; 
0.30 و‎ σου A items: 


and 


۷ 
= 
| V 
pa 


category C corresponds to 3 


V 
z 
V 
Ww 


category B corresponds to 5 


M 
{η 


312300 ۶۴ہ‎ 50۱۱45 το W > 


۳ 1۱۱۲ 1۸۵۲ 10۲5 OF THE TVC COST MINIMIZATION METHOD 

The UICP model assumes that the average annual variable 
costs are based on historical unit price and procurement lead 
time data and does not take into account that control over 
those two variables is possible during the procurement process. 
Price is, in fact, not a constant nor is lead time in a com- 
petitive bidding environment. In other words, the UICP model 
treats C and L as known values, and proceeds to solve for Q 
and R. However, when the procurement personnel solicit bids 
for reprocurement, Q is fixed and R has already been reached, 
while C and L are unknown until the vendors' bids are received. 

Additionally, the item manager is looking at only variable 
inventory costs, while the purchasing agent is primarily con- 
cerned with the purchase cost. The purchase contract is 
awarded to a "responsible contractor" who offers the item at 
a "fair and reasonable price" [Ref. 3] which generally means 
obtaining the material at the lowest bid price from among the 


vendors who can meet some required delivery date. 


270 


Reduction in the expected average annual total cost of 
stocking an item from that predicted by the "optimum" UICP 
solution is possible because negotiated variations in unit 
price and lead time can result in savings in some variable 
cost elements and/or the average annual purchase cost. A 
reduction in the latter usually outweighs increases in other 


variable cost elements. 


F. PARTIAL DERIVATIVES FROM THE TC EQUATION 

In order to understand their influence on the total cost 
equation, we will first take partial derivatives with respect 
to the lead time L and unit price C. The total cost equation 


is the sum of TVC and the average annual procurement costs; 





that is, 
TC = OC απο + BOC ۲ 
where 
E _ 4AD 
OC = ordering cost = a 
HC = holding cost = rc [Š + zo + B(O,R)] ; 
AE 
BOC = backorder cost = = B (O; R) ال‎ 
PC = procurement cost = 4DC 
Farst; 
(TC) 31000 , 3(HC) , 3(BOC) , 3(PC) 
91, oL aL, oL dL 
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and, since OC and PC are not functions of L, 








Bae remaining terms include B(Q,R), a complex function of L. 
The holding cost formula is a linear function of the 
safety stock zo and o is a function of L. The normal deviate 


1 not Q Function ofL since it is determined from the RISK 


formula. Therefore 
oH _ dO η. 0, Ε 
"πε" 


If we assume that we are considering Mark IV items and 


that MAD, is independent of L, then 


10857 MAD? -L EN MAD? - 


9g p) 1/2(1.57 MAD?) 


۱ 
ہے تشد 


2 
1.57 MAD, 


20 


Substitution of this formula into I gives 


2 
ee Z MAD, E 3B (Q, R) 


H _ 
a - τοι > οι τν 


Because B(O,R) is a very complex function of both C and L we 
will not attempt to obtain its partial derivative with respect 
to L analytically. Instead we will examine B(Q,R)'s behavior 


with respect to L empirically below. 


1 


Turning next to the BOC term, we get 


9(BOC) _ AE 3[B(Q,R)] 
dL 2 5 dL 


Substitution of these partials in the QTC equation gives 


IC(1.57 z MAD?) 


Ene D 3B (Q,R) AE 
0 20 + cn 
The first term of 2E is dependent on L in the denominator 


because of c and, since this is true, the term gets smaller as 
L increases. The behavior of this term is shown in Figure 1 


for the following set of parameter values. 


0 
lI 


8 units/quarter; 


E 


2.50 units/quarter; 


: 


QUartersS;‏ 222 و2 


3907 


E = 0.50; 
س0153‎ ٢ رت1‎ 


8.4-00: 


Its behavior appears almost linear in L and the rate of de- 
crease is quite small, being 0.07 on the average. 

The behavior of the second term of E can be deduced 
from the behavior of B(Q,R) as a function of L in FiGure n 
This figure shows that the relationship is virtually linear 


with a slope of 0.0012. 
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Thus, the second term is essentially a constant. The term 


9B (Q , R) 
Q 


the holding cost term (IC) and partially from the backorder 


mc + 5 which multiplies the term ۱۳ ۵ ۶۶ 
cost (45). 

Egure 3 suggests that the backorder term in the holding 
cost dominates the safety stock term for the values of L shown 
so that TC is essentially linear in L. In that figure the 
TC is called TVC since it does not include the fixed cost of 
annual procurement. 

The minimum value of L; namely, four quarters, shown in 
the figures corresponds to the current average administrative 
lead time value at SPCC. Since the procurement lead time 
includes both administrative and production lead time, realis- 
tic values of L can be expected to exceed this four quarters 
minimum value. 

Next we take the partial derivative of the TC formula with 
respect to the unit price C. The variables that change with 
C are Q, R, B(Q,R) and RISK. We can reasonably assume that 


A and o are not changing with C. 


AMI 3X0) 7 a (HC) 7 ð (BOC) 7 9 (PC) 
aC aC aC ος aC 
First we get 
ο 7 m 
>” G 4D 


EZE, the OC term contains Q which is a function of C at its 


optimal value. 


25 


26 


00 یں 


00 


0001 008 000 47 


el 


007 


0091 


0006 7 


0066€ 


ی 


۲ 10 10۵۲33۵1 e se 3500 A1O0JU9AUI SUT jo SqugSuOoduo;,) 


)5:171817110 (1 


üc 96 Fc GO 0c ul 91 Fi cI 01 0 
۱ ! ۱ ! ۱ ! ! ! ! ! 1 
O--29---O---O---AB--«- یبحم دی‎ ο πι ο - - T 


JA L V V 


JOH ο 


ου تک‎ 
)0 [1 





à e سم‎ = = u 
۱ = E سح‎ AY i : 
۰ را‎ 7 = = . 
۱ "eo 5 
B x سپ‎ “= k ۰ 
جو ھا ہے ه8‎ | 
E dA ce 





0t 


ποῦ. ο. 4λΡ.. . a [πο 
30 E 5c! Q ) = 4AD+ 56 8AD 


-1/2 I ADI 


—— = س ) 


l TC 
8AD 2 


= 4AD-5 (gap 


The holding cost term is a very complicated function of C since 
it depends on safety stock and B(Q,R). The safety stock is a 
function of z which is a function of RISK and hence a complex 
EIC On Of C. A plot of z as a function of C is shown in 
Figure 4 (the horizontal line up to a unit cost of 33 is a 
۰ {τοπ 1۱۶ 1 constraint that RISK > 0.01). The term 
B(Q,R) (which also appears in the backorder cost term of TC) 

is also a complex function of C and has been plotted in Figure 5. 


From these figures we might conclude that there is little 
OTC 
o 


of OC, HC, and BOC as a function of C in Figure 6 indicate an 


to be obtained from considering Fortunately, the plots 
almost linear relationship and we know that the PC term is 
indeed linear. Combining this information into a plot of TC 
as a function ofC, also shown in Figure 6, indicates that TC 
is virtually linear in C and is dominated by the PC term. 

in summary, we have considered the partial derivatives of 
TC with respect to L and C. We found that simple analytical 
formulas describing these derivatives could not be obtained, 
mainly because of the backorder term B(Q,R). However, empiri- 
cal analysis has shown that TC is virtually linear in both L 
and C. This fact will be important for the development of 


the technique to be presented in Chapter IV. 
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T Tp m SOC OST ANALYSIS. 


A. ANALYSIS METHODOLOGY 

In Chapter II we analyzed the behavior of the expected 
۱39 ۰ annual costs (TC) as a function of C and L by taking 
the partial derivatives with respect to each and studying 
their behavior. This provided an indication of the indepen- 
dent influences of C and L on TC. That analysis and the 
analyses to be presented in this chapter are based on the 
important assumption that TC represents a steady state situa- 
tion in which the bid C and L values being evaluated will be 
the same for all future reprocurements. In this chapter we will 
ENEN zc the behavior of TC when C and L are both allowed to 
vary. In particular, we will attempt to determine those C and 
L combinations which give the same value of TC. 

Changes in both C and L from those currently in the ICP 
computer files are expected from a reprocurement buy. In fact, 
we expect each vendor who bids on the procurement contract will 
Submit a bid having C and L values which are different from 
those of the other vendors. 

In Chapter II we realized that the partial derivatives were 
complex functions of C and L and that we had to resort to 
computer-aided analysis of a simple numerical example. It is 
clear that allowing both C and L to vary will make the analysis 
problem even more complex and therefore we will also use 


computer-aided analysis of that same example in this chapter. 


رد 


The computer program for this analysis is contained in 
Appendix A. It is designed to calculate, for each combination | 
of C and L, the order quantity (Q), the reorder level (R), the 
expected number of backorders and then the values of the 
components of TC. All of the formulas and constraints from 
Chapter II are included. As in Chapter II, MAD. is assumed to 


be independent of the bid value of lead time, permitting the 


same values to be used in computing o for different bids. 


B. VALUES OF C AND L 

Both the ranges of C and L and the increment within each 
range need to be selected before the analysis can beicon DEE 

The range Of lead time values we will use is from four to 
16 quarters. The four quarters lower bound corresponds to 
the current average value of the time between when the reorder 
point triggers a buy and the contract is signed with the vendor. 
The maximum lead time values we have observed have been about 
4 years or 16 quarters. The lead time increment used in the 
computer program was selected to be 0.25 quarters as a compro- 
mise (daily increments were considered to be much detail and 
quarterly increments were considered to provide insufficient 
detail). 

The selection of a range for C is complicated by its role 
with the expected quarterly demand, D, in defining when cer- 
tain other UICP model parameters change. They are used in 
defining the borders between Mark II and Mark IV and in defining 


the border between different values of the administrative 
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order cost A. Crossing the border between Mark II and Mark 
IV requires changing the equation for σ. nat ۱۹71 18 406 
by CD = 75. For our example, we assumed the quarterly demand 
to be 8 so the border between the marks will be at C - = = $9.375. 
When the total purchase cost of a buy (CQ) exceeds a cer- 
tain volume the value of A changes. At the time of the initial 
work done on this vendor bid problem [Ref. 5] the value of A 
changed from $380 to $1050 when CQ exceeded $8000. These 
values have also been used in our example. The 1985 values 
have increased to $390, 1080, and 10,000, respectively. 
In our example we don't want A to change from the assumed 
value of $380, therefore we need to find the range of C which 


ensures this for D = 8. 


Now 
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lun r rore, if CO « 8000, then 


325/C « 8000 


resulting in 


Based on these "boundary" values for C, we selected a 
range for C, from Sl up to S600, and we will use increments 


er Sl. 


C. TOTAL COST CURVES 

The computer program was designed to run the ranges of. 
unit costs from $1 to $600. For each unit cost the program 
computes a TC value for each lead time value from 4 quarters 
to 16 quarters in steps of 1/4 quarter. A sample printout is 
shown in Appendix B. However it is easier to visualize and 
study these results if they are presented graphically. 

Ideally the results should be shown in three dimensional 
graphs. Unfortunately it is very difficult to present and 
understand the results from that kind of graph. Instead we 
have chosen to show the results on two dimensional graphs. 
Figures 7 and 8 show the results. Every 48 points on the x 
axis correspond to one unit cost value and all of the 48 lead 
time values ranging from 4 quarters to 15.75 quarters. This 
axis is called 'C and L' to reflect the sequencing. For exam- 
ple the first graph (Figure 7) shows a series of sloping steps; 
the first at the origin has a unit cost of Sl and the x axis 
values from 0 to 47 correspond to all the L values. The TC 
values of this first step begins $181 for L = 4 quarters and 
gradually increase as L increases to a value of $197 for 
L = 15.75 quarters. The second step corresponds to C = $2 
and its TC values range from $233 0م‎ $258. 

Figure 7 shows how TC behaves for very small values of C. 
From $1 to $9 the item is a Mark II item. Above $9 it is a 
Mark IV. The change from Mark II to Mark IV occurs at the 
"C and L" values of 432. Figure 7 has a maximum σ valucmem 


$30. Figure 8 shows C values from $21 to $40. 
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The behavior of the TC curves in Figures 7 and 8 for the 
Mark II and Mark IV ranges of C are similar. For the low C 
values in each range, each TC step is higher than that for the 
preceding C values for all of its range of L values. Thus, if 
we had two bidders whose C values are in this range, the bidder 
with the lowest C value would be the winner regardless of his 
L value. 

On the other hand, after C passes a certain threshold for 
each Mark, we see the same TC values can occur for two adjacent 
values of C. In fact, Figure 9 shows this can occur for more 
than two C values. In such a situation, we would not want to 
automatically select a bidder with the lower C value until we 
had compared his L value with that of the higher cost bidder. 
The threshold C value for Mark II appears to be $5 and, for 
Mark IV, it appears to be $38. 

The reason for the change in slopes between the C = 9 curve 
(Mark II) and the C = 10 curve (Mark IV) is due to the differ- 
ence in the formula for the standard deviation of lead time 
demand between the two Mark codes. Figure ll shows that the 
Mark II standard deviation increases with L much more rapidly 
than for the Mark IV. Since R and hence B(Q,R) are increasing 
functions of L because of their relation to the standard 
deviation, the shortage costs and safety stock holding costs 
terms of TC increase more rapidly with L for Mark II than for 


Mark IV for any given C value. 
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Pee ISOCOST LINES 

Isocost is a term from microeconomics which means that we 
have numerous combinations of two resources which can be used to 
produce the same quantity of product. As Heinz Kohler observed 
[Ref. 4], "We can draw a family of straight isocost lines, each 
of which shows all the alternative combinations of two inputs 
that the firm is able to buy in a given period at current 
market prices, while incurring the same total cost." As we can 
see from Figures 7, 8, 9, and 10, the variables of inter- 
est, C and L, do have several combinations which produce the 
same TC value. The questions which remain to be answered are: 


l. Can isocost curves be constructed for all L and C 
values in the ranges we are considering? 


2. What is the shape (linear or non-linear) of those 
isocost curves that can be constructed? 


The answer to the first question is that isocost curves can 
only be constructed for those C values above the thresholds 
observed in the last section (C = $5 for Mark II and C = $38 
for Mark IV). The answer to the second question requires us 
to examine some C and L pairs. 

Figure 12 presents thirteen combinations of C and L in 
the Mark IV category. These combinations correspond to a 
nominal value of TC of $13,500. In reality, searching for 
C, L pairs which give TC = $13,500 precisely is extremely time- 
consuming. What is shown are C,L pairs which gave TC values 
within the range of $13,500 + $5. A least squares fit of a 


straight line through these points resulted in a correlation 
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coefficient of 0.999 and therefore the shape of the isocost 


née can be considered as linear. The equation for the isocost 


line is 

L = -1.71C + 520.29 
Figure 13 shows the resulting line. This figure also shows 
the line for TC = $13,510 and indicates that as TC increases, 


the isocost line would move to the right. 

A computer program (Appendix C) was written to compute the 
parameters of the isocost lines for TC values from $6800 to 
$13,400. For each TC value C,L pairs were used if they re- 
sulted in a value which was within $5 of the desired TC value. 
Typical results from the program are presented in Appendix D. 
From these results a curve of the slope a was plotted as shown 
in Figure 14. If the oscillations are ignored, it is clear 
that the slope becomes less negative as TC increases and that 
theoretically a family of isocost lines should appear as shown 
in Figure 15. 

The oscillations are a consequence of the use of C,L 
combinations which are not exactly the value of TC. In other 
words, the inaccuracies induced by "rounding" to the nominal 
value of TC create least square fits with slopes which oscil- 
late. Table 1 illustrates the situation. As the rounding 
range is reduced, the oscillations are also reduced and, in 
the limit, would provide a smooth curve of a as a function 


EL TC. 
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Figure 15, Theoretical Shape of Isocost Lines 


We could fit a curve to the data of Figure 14 to obtain 
the function a = £(TC). However, we are more interested in 
providing a tool for procurement personnel and therefore we 
prefer to discuss a curve of the slope a as a function of C. 
Each TC value has a certain corresponding value of the slope a 
as shown in Figure 14. On the other hand, each TC has a range 
of C which creates this TC as shown in our example in Table l 
and theoretically in Figure 15. We can compute the average of 
these C values from the range of values used in the computer 
program to generate Figure 14 and relate this average to the 
corresponding slope values via the corresponding TC value. 

The graph of the slope a versus C which results is shown in 
Figure 16. As should be expected, this curve also has 


Oscillations. 
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We can now attempt to fit a mathematical curve through the 
points on this graph. Although we could try a variety of. 


formulas, our intuition suggests that the curve behaves as: 


l 


aC +B‏ ات 


We see that the denominator is a linear function of C. The 
parameters o and 8 can be found by using two extreme points 
of the curve. For this example, when C = $178 the slope a 
is -2.72, and when C = $302, the slope a is -1.7. Substituting 


these values to the above formula, we get 


1 


4.14 = 178 +B 


lo 
3020 +8 


which can be rewritten in the following form 


484.16% + 2.728 = -l ; 


SSO IG = -l 


Solving these two equations for the two variables gives: 


0 0536625” 


e 
II 


= 0 0017655 


[ov 
II 
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The formula for the curve will therefore be: 


1 
-0 .001 7/638 2 00536 am 


This curve has been plotted on Figure 16 and we see that it pro- 
vides a very good fit to the points on the graph. While more 
precise methods can be applied to find the curve, we suspect 
that in most cases the function we have selected reflects the 
general shape of the curve. 

From this curve we see that for small values of the unit 
cost the slope a is very negative which means that a large 
reduction in L is needed to compensate for a small increase 
in C in order to maintain the same total cost value. For 
large values of C the slope a approaches -1 which means the 
isocost line has an angle of -45°. Thus, the unit cost and 
the lead time have almost the same influence. A small increase 
in the unit cost can be compensated for by a small reduction 
in the lead time in order to remain at the same total cost 


value. 
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IV. A BID EVALUATION TECHNIQUE 


A. INTRODUCTION 

It is obvious that anyone could use a computer program 
similar to that used in Chapter III to compute the TC for each 
bid. However, this would require that a procurement officer 
have access to the ICP mainframe computer so that he could make 
several computer runs. This could take considerable time and 
money even if computer capacity was available. Because this 
is not usually possible at an ICP due to saturated computer 
capacity, some simpler approach is needed for bid evaluation. 
The purpose of this chapter is to present such an approach for 
consumable items. This approach would require, at the most, 
access to a micro computer. In most cases, it would require 


Pisa hand-held calculator. 


B. BID EVALUATION 

The bid evaluation procedure is shown as a flow diagram 
in Figure 17 and will be discussed in this section. It can be 
Meea τος any number of bids but can consider only two at a time. 
In the following discussion we will therefore consider only 


two bids, denoted by (C4 , L4) and (C L5). 


2 


The first step is to determine which bid has the lower unit 
cost C since this is the dominant bid parameter. Three cases 


Ean Occur: C1 = Cos Ci < Cor Or Ci > C: 
If we have Ci = C. then we must examine the values of Lj 


and Πο. T£ L; = L, then the two bids are identical in the 
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Figure 17. Bid Evaluation Flow Diagram 
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sense of the bid parameters and vendor selection is left to 


the discretion of the procurement personnel. If L; < La y then . 


(61,1) is the better bid since it provides the lower TC 
value (as Figure 7 showed). If L, < Li, then (C5, L5) is the 
winner for the same reason. The upper right side of Figure 17 


describes these steps. 


Hr Ci < C, and if Li 


better. Similarly, if Ci 


E bo ΙΕ; 122 (C4 , I4) is the 


2 


> C, and L, < L4 


the better. Again, the reasoning can be confirmed by referring 


then (C5, D) ls 


to Figure 7. The steps for these comparisons are shown in 
the top part of the left side of Figure 17. 
1 2 and if Ly > L, or C. < Ci and L, > Lı then we 


must determine which bid has the lower TC value by using o and 


ECS <C 


6. We first compute the average of Ci and Cai that is 


We then compute the slope associated with C using the follow- 
Fi formula: 
1 


a - [ας  β] 


Next we compute the differences in the C and L values using 


the following formulas: 


ENE C 
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AL. = OL. 


Finally, we compare AL with the value of the product a*AC. 


If AL > a*AC then (Ci,L is the better bid; if AL ES 


1) 


then (C.L is the better bid, and if AL = a*AC both bids 


2) 
are equally good. The steps of this comparison are shown in 
the bottom part of Figures” 


The reasoning behind these comparison steps can be deduced 


by considering Figures 18 and 19. 





Figure 18. (Cj,,L 1) with Figure 19 (C1,L4) with 
Lower TC Value (1 017 اه‎ τὸ Valu 
In Figures 18 and 19, we see the isocost line for C passing 
through the point (C,10) where L = 10 is the average value of 
L. The value of a* is, in fact, the slope of the isocosSte m 
passing through this "average point" as a consequence of the 


analysis from the last section ot Chapters 
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In Figure 18 we see that (Ci , I4) has a lower TC value than 
(C5,L5). Em cmcomecEEnesectwo"points"with the dotted line, 


the slope of that line is: 


۵ 

| 
D>] > 
ale 


and we see that a » a* since it is less negative. Now we can 


rewrite that inequality as 


and, after multiplying both sides by AC, we obtain the first 


result that 


ἂν و‎ 


In Figure 19 the (62,12) bid has the lower TC value. In this 


case the slope of the line is 


۲ 271. A 
Ci - C^ ΕΙ AG 
In this case, however, a « a* since it is more negative. Thus 


Nb η " 
AC = a < a 


and therefore 


AC‏ ص ان 
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C. AN EXAMPLE 


The example in this section illustrates the technique des- 


cribed above. Let us suppose that the current price and lead 
time are $230 and 10 quarters, respectively, a = ο ου το un 
B = -0.0537, and the rest of the parameters are those of the 
example in Chapter II. Now suppose we have received four new 
bids: 

a) $228 and 15 quarters; 

b) $232 and 5 quarters; 

c) $232 and 6 quarters; 

d) $233 and 4 quarters. 

As we explained in Section B, we will consider two bids 

at a time. Let us take b and c. In this case we denote them 


as l and 2, respectively. 


5232 PA SR 


Ω 
I 


$232 L, = 6 


© 
| 


Figure 17 tells us that when we have σι = Cor we have to check 


L4 and L, - In this case L, 


better. Now let us compare bid d with the winner bid (b). 


is greater than L4 so bid 1 is 


As above, we now denote these as l and 2, respectively. 


S293 He M c 


Q 
II 


$232 T. EE 


(2 
I 
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Now we must use the left side of Figure 18 because we have 


C, is greater than C. and L, is greater than L We have 


1 1 


to compute the slope a*. The first step is to compute C. 


C = a = 232.5 


1 ۱ ۱۱۳20 a, 8, and C values into the a* formula, we obtain 


Ce ORO 1 ۰ = 7 
Next we need the value of AC and AL. 
AC = C-C = 232 - 233 = -1 
AL = L -L = 5- 4 = 1 
The AL test needs a*AC. 
a*AC = (-2.16)(-1) = 2.16 


comparing AL with a*AC, we see that 


nr AS, 


and therefore (C is the better bid. This corresponds to 


2132) 
pud b. 
Now we must compare bid a to bid b. We denote them as 1 


and 2, respectively. 


0 


$228 L = | See 


Ω 
| 


5232 n = گے‎ 


O 
I 


In this case again we must use the left side of Figure 17 


since C. is greater than C, and L is greater than L 


1 2" 
چ‎ = 2284232 . jq 
a* = ([-0.00176-230 « (20.053)]72 . 
AC = C3 = Ci = 232 = 2283 = a 
AL * LQ,- Lh, = 5-15 = -10 
a*AC تم 4( 2 کان‎ « <7 


Comparing AL with a*AC we again have 


AL στο 


and therefore (C5, L5) is better which means that bid b is 
a better bid than a. We also know now that it is the best of 
the four bids. 

We can now compare this result to the real total costs 
from the computer program as shown in Table 2. Table 2 shows 
the TC values for each bid determined by the computer program 


and- -confirms GUI Conclusion: 
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TABLE 2 


Computed Values of TC 


Bid Unit Price Lead Time Total Cost 
G L 
a 228 ڑا‎ 10,700 
b 232 5 Io De best) 
E 232 6 LO 700 
d 233 4 10,700 
current terms 238 10 10,690 


D. MODIFICATION TO THE UICP PROGRAM 

۳ ۵11 1 ۵۱۳6 this approach two additional item parameters 
need to be computed and stored in the Master Data File (MDF). 
These should logically be computed each quarter by the cyclic 
levels and forecasting application (D 1). The additional 
parameters are the values of oa and B from the curve fit to the 
plot of a as a function of C (see Figure 16). The subroutines 
needed to make these calculations can be developed from the 
programs in Appendices A and C. Then, when the request to make 
a buy of an item is sent to the procurement personnel by the 
inventory manager, the former can use these parameters in 


evaluating bids from vendors. 
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V. SUMMARY AND CONCLUSIONS 


A. SUMMARY 

Chapter II presented a brief overview of the purpose and 
underlying assumptions of SPCC's UICP wholesale consumable 
procurement model. The model's total variable cost equation 
and optimization results were presented, and it was shown that 
the model determines the optimum order quantity and reorder 
level for an item based on that item's forecasted quarterly 
demand rate, procurement lead time and unit price. Later in 
the chapter we analyzed the behavior of the total inventory 
cost as independent functions of the unit price and the pro- 
curement lead time. The total cost was found to be a linear 
function of each of these two variables. 

Chapter III examined the behavior of the total costs as a 
combined function of these two variables. The results showed 
that equal values of the total inventory costs existed for 
different pairs of C and L. This fact suggested that we could 
find an isocost curve describing all the combinations of C 
and L. Subsequent analysis showed this curve to be a straight 
line. The slope of this line was also found to be a 
nonlinear function of unit costs and a formula was found which 
described that یز‎ 

Chapter IV then used this information to develop a simple 


technique for evaluating different bids. 
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B. CONCLUSIONS AND RECOMMENDATIONS 

ocd een vent@nyscOceS can be reduced through modifications 
of the procurement process to include procurement lead time as 
well as unit price in selecting the winning bid for stock 
replenishment contracts. The simple technique developed in 
Chapter IV can be used to efficiently evaluate bids containing 
lead time as well as unit price. 

Implementation of the technique can be done in three major 
steps. 

The first step would be to add a subroutine to the UICP 
"levels" application (D 1) which would calculate for each item 
the parameters needed for using the technique. After calcu- 
lating these parameters (which should be done quarterly), they 


should be included in the Master Data File for each item so 


that they can be given to procurement personnel when a replenish- 


Heme buy is needed. 

The second step would be to develop a procedure for provid- 
ing the procurement personnel with the additional information 
and computing equipment needed to use the technique. 

The third step would be to train procurement personnel in 
the use of the technique. It may also be appropriate to ex- 
plain the evaluation technique to potential vendors who are 
going to compete for supplying the item to the Navy. Some 
initial monitoring of the use of the technique should also 


be considered. 
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C. AREAS FOR FURTHER RESEARCH 

The fit of a mathematical expression to the curve of the 
slope of isocost lines as a function of unit price needs 
further investigation. Statistical validation of the proposed 
mathematical expression from Chapter III is needed and possible 
other functions should be investigated also. 

In the implementation of the additions to the UICP levels 
application D 1, discussed in Chapter IV, several modifications 
of the programs from Appendices A and C should be investigated 
which have the potential for reducing the computational steps. 
One modification woulá be to limit the range of unit price 
values used to determine the parameters needed for the bidding 
technique to those values which would be expected to be bid. 
Perhaps, for example, a range of unit price values from ninety 
to two hundred percent of the current value stored in the Master 
Data File is reasonable. 

Another related improvement would be a procedure for auto- 
matically selecting the incremental value of unit price to be 
used to develop the curve of isocost slopes as a function of 
unit price. The size of the increments influences the computer 
program's ability to compute isocost slope values and the time 
required to determine them. Manifestations of the problem are 
the "threshold values" observed in Figures 7 and 8 of Chapter 
III. Below such values it appeared that no slopes could be 
or needed to be computed. In reality, these threshold values 


go to zero as the increments of the unit price go to zero 
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and slopes can indeed be determined for all positive unit 


prices. However, as a unit price approaches zero, the iso- 


cost slope approaches negative infinity. A tradeoff analysis 


needed between the size of the increment and 
isocost slope values for small unit price 
results of such an analysis should provide a 


"rule of thumb" for determining the unit cost 


1s therefore 
the need for 
values. The 


formula or a 


increment size to use for developing the isocost slope 
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Curve. 
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۱۳۲ IS A NONINTERACTIVE PROGRAM UTILIZING THE WAIFIV COMPILER WHICH * 


INPUTS AN INVENTORY ITEM'S CHARACTERISTICS FROM A DATA FILE AND THEN 
DETERMINES THE OPTIMAL ORDER QUANTITY AND REORDER LEVEL AND THE 
AVERAGE ANNUAL COST OF STOCKING THE ITEM . THE PROGRAM GOES OVER 
NUMEROUS VALUES OF UNIT PRICE (C) AND LEAD TIME (L),AND FOR EACH 
COMBINATION CLCULATES THE ABOVE PARAMETERS.THE PROGRAM IS ACTIVATED 
BY AN EXEC FILE WHICH DEFINES THE OUTPUT FILES AND LOADS THE IMSL 
SUBROUTINES. 


e de de de de de de e de de Gc e de ce de de ce de e e de Y jc jc $c de de de cde e de ce de ας de de ce e de de ce ας de de de cde ce ce cde de cde de cde cde de de de de de de ce de e de de Tc $e eoe‏ سس 


te te te te ας de de de ας te VARIABLE DEFINITIONS de de de te te te te He te kx 

A = ADMINISTRATIVE COST OF PLACING AN ORDER 

BOC = AVERAGE ANNUAL COST OF BACKORDERS 

G = UNIT PRICE 

ECC = INITIAL UNIT PRICE 

CQ = PURCHASE COST 

AVERAGE QUARTERLY DEMAND RATE 

ITEM ESSENTIALITY (MILITARY WORTH) 

= EXPECTED NUMBER OF BACKORDERS JUST BEFORE AN ORDER ARRIVES 
= ERROR NUMBER 

REQUISITION FREQUENCY (D/S) 

MAXIMUM BOUND ON Q 

AVERAGE ANNUAL HOLDING COST 

= INVENTORY HOLDING COST RATE (FRACTION OF UNIT PRICE PER YEAR) 
COUNTER FOR UNIT PRICE 

COUNTER FOR LEAD TIME ۱ 

SOURCE OF INPUT DATA-DETERMINED IN VARIABLE DECLARATION 
COUNTER FOR 'C AND L' 

L = PROCUREMENT LEAD TIME 

LAMBDA = STOCKOUT COST RATE ($/UNIT/YEAR) 

MADD = MEAN ABSOLUTE DEVIATION OF DEMAND 

MADL = MEAN ABSOLUTE DEVIATION OF LEAD TIME 

OC = AVERAGE ANNUAL ADMINISTRATIVE ORDERING COST 

OUT = DESTINATION OF OUTPUT. DATA=DETERMINED IN VARIABLE DECLARATION 
P = UNCONSTRAINED RISK 

PC = AVERAGE ANNUAL PURCHASE COST OF THE ITEM 


B 
R 


II II! XA O Il M 


E 


D 
2 
E 
E 
El 
G 
H 
I 
I 
ET 
IN 
K 


۱ ۱ ۱ Il 


PMAX = MAXIMUM RISK CONSTRAINT 
PMIN = MINIMUM RISK CONSTRAINT 
POUT = CONSTRAINED RISK 


PPV = PROCUREMENT PROBLEM VARIANCE 

Q = UNCONSTRAINED ORDER QUANTITY 

QHAT = FINAL CONSTRAINED ORDER QUANTITY 

RHAT = CONSTRAINED REORDER LEVEL 

ENECAVERAGE NUMBER OF UNITS PER REQUISITION 

SIGC - STANDARD DEVIATION OF MEAN LEAD TIME DEMAND 

TC = AVERAGE ANNUAL TOTAL COST OF STOCKING THE ITEM 

TVC = AVERAGE ANNUAL TOTAL VARIABLE COST 

4 = PROCUREMENT PROBLEM VARIABLE (MEAN LEAD TIME DEMAND QUANTITY) 


*** SENTINEL VALUES FOR GENERATING ERROR/WARNING MESSAGES TO THE USER 


ERR, WAFNCO, WARNO, WARNR 
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ΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩΩ 


E 


E * κκ κ e چا‎ e eode * * T Xx VARIABLE DECLARATIONS * e e e e oe e oe ΧὩ de oe ve κ κ hee deo ox 


C 
REAL A,C,D,E,F,H,I,L,LAMBDA, MADD, MADL,OC,P,PC,PMAX, PMIN EOS 
*PPV,0,QHAT ΑΕΙ ΑΡ S 16 ٢ 1.7 
INTEGER ERR, IN, J,OUT, WARNCQ,WARNO, WARNR 
DOUBLE ۰ 7+ B, BOC EBO 
E 


DATA IN/4/,0UT/6/ 
10 CONTINUE 
C *** READ INPUT DATA FROM FILE *** 
20 READ (IN,450) A,C,D,E,H, I, L, MADD, MADL,S 
C *** CHECK FOR END-OF-FILE SENTINEL VALUE IN FIRST COLUMN OF INPUT*** 


2 IF (¿NOT A GE 100) ο τοι... 
C *** WRITE INPUT DATA TO TERMINAL *** 
ccc=c 
KQ=11040 
SL=L 


C *** DO LOOP OVER VALUES OF UNIT COST 
DO 902 112281 250 
GST 
C *** CACULATING ORDER QUANTITY WITH CONSRAINTS 
QSSQRT(8*A*D/(I*C)) 
CAMA TDI) 
OCON=12*D 
OHAT=AMIN1(G,QCON) 
C *** DETERMINING ADMINISTRATIVE COST 
CQ-C*QHAT 
IF(C*QHAT.LE.8000.)GO To 54 
A=1050.00 
54 CONTINUE 
C *** DO LOOP OVER LEAD TIME VALUES 
DO 901 II=1,48 
SL=4+(11-1)/4. 


KQ-KQ*1 
a 
C *** DETERMINE ITEM COG/ASSOCIATED VALUES,CALCULATE PROCUREMENT PROBLEM 
C *** VARIABLE AND PROCUREMENT PROBLEM VARIANCE *** 
C : 
30 CALL DATACK (C,D,F,1I,SL,AMBDA, MADD, MADL, PMAX, PMIN, PPV, S,2Z, ERR, 

* OUT) 


C *** GENERATE ERROR MESSAGE IF DEMAND DATA DOES NOT MEET REQUIREMENTS 
C *** FOR ASSUMPTION OF NORMAL DISTRIBUTION OF LEAD TIME DEMAND QUANTITY 
GO TON Sd 82 39 1) ۷۳ 


e WRITE (OUT,600) 
co πο lli 
32 WRITE (OUT,610) 
CONTO A 
33 WRITE (OUT,620) 
GO TO II 
34 CONTINUE 
e 
C *** CALCULATE/CONSTRAIN RISK *** 
e ۱ 
40 P - (D*I*C) / ((D*I*C)*(LAMBDA*F*E)) 


IF (.NOT.P.GT.PMAX) GO TO 41 
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POUT = PMAX 
«ο το 43 
IF (.NOT.P.LT.PMIN) GO TO 42 
POUT = PMIN 
CO 1O 43 
POUT = P 
CONTINUE 
SIGC=SQRT( PPV) 
CDFA1=1.-POUT 
CALL MDNRIS(CDFA1,2ZNOR, IER) 


NOTES ٧11815 IS THE CEST OMEUTER CENTIER IMSL ROUTINE FOR THE 
INVERSE NORMAL PROBABILITY DISTRIBUTION 


RHAT=ZNOR*SIGC+Z 
CALL EBOCAL(EBO, QHAT, SIGC, RHAT, Z, BRRQ) 
B=EBO 


COMPUTE THE AVERAGE ANNUAL TOTAL VARIABLE COST, AVERAGE *** 
ANNUAL PURCHASE COST, AND AVERAGE ANNUAL TOTAL COST κ κκ 
ΕΜΕΟ ΟΕ ΕΕ ΕΕ کلانا لصا ناد رل‎ ΟΠΑΤ ΚΕπατ τνς οὐ πς, σος, 5) 


PC-C*4*D 

TC=TVC+PC 

TCI=TC/10.+0.5 

ICII=AINT(TCI)*10. 

ICIZ=AINT(TC) ۰ 

IC=C 

1۳ ٢ 

WRITE (OUT,750)IC,SL, ICII,QHAT,RHAT,KQ,A 

FORMAT(2X,I4,1X,F5.2,2X,I5,2X,F5.1,2X,F6.1,2X, I5,2X, F5.0) 

CONTINUE 

CONTINUE 

CONTINUE 

ΤΟ 0 

CONTINUE 

STOP 

ο n τον EG O JX 77.2 41X/F4.2,1X,F4.1,1X,F4.2,1X, 
BS e εδ 2. 1X. ος ο 1X,E5.2) 


FORMAT ('O','INVALID DATA - D IS LESS THAN 0.25 = LEADTIME 


*DEMAND NOT NORMALLY DISTRIBUTED') 


FORMAT ('O','INVALID DATA - F « 3 AND Z < 4 = LEADTIME DEMAND 


*NOT NORMALLY DISTRIBUTED') 


FORMAT ('O','INVALID DATA- F « 1 AND Z « 20 - LEADTIME DEMAND 


*NOT NORMALLY DISTRIBUTED') 


END 


41 


42 


C *** 
C «κκ 


CQ «κκ 
C kkk 
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gol 

902 
(11 
120 
T30 


450 


600 


610 


620 


€ 


C*** END OF MAIN PROGRAM AND BEGINNING OF SUBROUTINES *** 


SUBROUTINE DATACK (C,D,F,I,L, LAMBDA, MADD, MADL, PMAX, PMIN, PPV,S,Z, 
ERR, OUT) 


e 


CC * Z * x * k k k k k k k k k k k k k k k k k k k k k Kk k k k k پو پو پو پو‎ k kok k پو پو پو پو‎ KOR پو پو پو پو پو‎ eoo koe ck oko x kk oko ko ko ox xo X 


2 THIS SUBROUTINE DETERMINES THE ITEM'S COG, AND RETURNS THE SHORTAGE- 
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COST FACTOR, MAXIMUM ALLOWABLE RISK OF STOCKOUT, THE PROCUREMENT 5 
PROBLEM VARIABLE AND THE PROCUREMENT PROBLEM VARIANCE TO THE MATE 
PROGRAM UICP1. IT ALSO GENERATES AN ERROR MESSAGE IF THE LEADTIME * 
DEMAND DISTRIBUTION IS NOT NORMAL (BASED ON THE ۲۱ ال‎ M * 


تر ہی ری ہی ہی ہہ ہی عا غاب چب غاب عب اب عاب چب Ñ W W‏ کا Wc 7 W W W Wc T W T‏ غاب KKK KKK ERR‏ غاب غاب عا غاب غاب وی و و و و وس W W W W W T T‏ ہہ یں 


VARIABLE DECLARATION *** 

INTEGER  ERR,OUT 

REAL  C,D,F,I,L,LAMBDA,MADD, MADL, PMAX, PMIN, PPV, S, VS, Z, SIG2D, SIG2L 
A 

⁄ s 


Mr N 
Ail اا‎ 
e Ü O 


R 


IF ( NOT DILE "۰ ۶+ 7٤( 


ITEM IS MARE CODE GOS 
ERR Z1 
GOTTO 350 


TF (Nol 29 1 Se CO TO 240 
ITEM IS MARK CODE 1 OR 3 *** 
PPV = (2.028 = 2 2ک 701ر‎ 
CO TO 220 
CONTINUE 
TE NOT eD LE SI CO TOO 
ITEM IS MARK CODE 2 *** 
SIG2D= 1.57 * MADD * MADD 
SIGE L 
PPV = L*SIC2D + D*D*SIG2L 
GO TO 220 
CONTINUE 
ITEM IS MARK CODE 4 ۶ 
51629 = 1.57 * MADD * MADD 
SIG2L = 1.57 * MADL * MADL 
PPV = L*SIG2D + D*D*SIG2L 
CONTINUE 
EXCESSIVE VARIANCE SCREEN *** 
VS = PPV / Z 
IF (-NOD.VS.CT.150, CO το ο 
ΕΞ 2. 112 * 271 1 R 
CONTINUE 
CHECK DISTRIBUTION OF Z, ASSIGN LAMBDA,PMAX,PMIN BY ITEM COG *** 
IF ενοτσο  }οο το 970 
ITEM IS CATEGORY A *** 
LAMBDA = 1500.00 
PMAX = 0.30 
PMIN = 0.01 
GO TO 340 
IE (NOTE ESE Na.) G0 0 
ITEM IS CATEGORY A *** 
LAMBDA = 1500.00 
PMAX = 0.30 
PMIN = 0.01 
οσο πο ٤ 
IPC- NOD E CEI econ Tomaso 
ITEM IS CATEGORY B *** 


0000000 


κ‏ ده دا 


C xxx 


100 


C xxx 


200 


C **x 


ن 21 


کب ځا E‏ 


220 


Qe 


230 


#ښ ىے٭ C‏ 


300 


C «κκ 


310 


ο x x 


320 


E kkk 


LAMBDA - 1000.00 
PMAX = 0.40 

PMIN = 0.01 

IE (Z.LT.4.) ERR = 2 
GO TO 340 


330 CONTINUE | 
۰ 6 *** ITEM IS CATEGORY C *** 
LAMBDA = 500.00 

PMAX = 0.50 
PMIN = 0.01 
τε Le eco.) ERR = 3 
340 CONTINUE 
350 CONTINUE 
RETURN 
END 
SEJECT 
C 
SUBROUTINE TVCOST (A,B,C,D,E,I,L,LAMBDA,QHAT,RHAT,TVC,OC,HC,BOC,S) 


په چب اب چک < چې چا اب * اب اب ار اب * اب اب اب اب * * * — * * de cde cde‏ — * < — * < * اب ede e cde de cede dece cede * K‏ اګ (C * * * * * * * * < RARA KR e‏ 


THIS SUBROUTINE CALCULATES THE ANNUAL TOTAL VARIABLE COST OF 5 
ο. Οκ THE ITEM PER THE UICP CONSUMABLES INVENTORY MODEL, τ 
. HSIBERING TIME-WEIGHIED,ESSENTIALIIY-WEIGHIED REQUISITIONS 2 
۰۰۰۷ ΙΝ ACCORDANCE WITH DODINST 4140.42. = 
QC ke de پر پر پر پر پر پر پر پر ې ری‎ Fe پا پا پا پا پا پا پا پا پا پر‎ eoe kc پا پا پو‎ e پو پر پو پر پو پا پو پر پو پر پر پر‎ RRR RR RRR RR KR KE 
C 

ENS VARIABLE DECLARATION *** 

ο Σο PRECISION BOC BEBO 


REAL A,C,D,E,I,L,LAMBDA, QHAT, RHAT, TVC,OC, S 
OC = A*4*D/QHAT 
HC = I*C*(QHAT*.5+RHAT-D*L+B) 
BOC = B*E*LAMBDA/S 
mic OC + HC + BOC 
RETURN 
E 
END 
SEJECT 
SUBROUTINE EBOCAL (EBO,QHAT,SIGC,X, Z, BRRQ) 
Οκ κκ. e de de de de eoe ode dede de de de eee ode cde ode oe اب اب‎ de de de cde dede oe oe dede cde de RT KKK RRR RRR RRR RRR چاو‎ RRR KKK RK KK 
5 THIS SUBROUTINE CALCULATES THE EXPECTED NUMBER OF BACKORDERS AT 
C THE END OF THE CYCLE FOR A PROPOSED ORDER LEVEL (X), GIVEN A LEAD 
C TIME DEMAND WITH MEAN OF Z AND STANDARD DEVIATION OF SIGC. 
C * * * * k < * * چا — ې ځا اب اګ * )ې * * * * * اب اب * * اب‎ RETR o o ve ال چل پل چې‎ e ve پل‎ k k k T k — ځا * < * — < * < * — چب پب‎ —< — = 
C *** NOTE - MDNORD IS THE NPS COMPUTER CENTER IMSL ROUTINE FOR THE ۶ 
C *** NORMAL PROBABILITY DISTRIBUTION OF A DOUBLE PRECISION VARIABLE *** 
C 
C *** VARIABLE DECLARATION *** 
DOUBLE PRECISION  A1,Bl1,CDFAl,CDFBl1,CCDFAl,CCDFBl,PA1,PBl,PI,EBO 
REAL X,SIGC,Z,QHAT 
= 3.1415926535 
AZ) 756 
= (X + QHAT - 2) / SIGC 
PAl-1./DSQRT(2.*PI)*(DEXP(-A1*A1/2.)) 
PBl - 1./DSQRT(2.*PI)*(DEXP(-B1*B1/2.)) 
CALL MDNORD (Al,CDFA1) 
ECLEAÎ = 1 - CDEA1 


68 


CALL MDNORD (B1,CDEB1) 

CCDEB1 = 1. - 1 

BR=( (SICC*SIGC+(X-2Z)*(X-2))*CCDEA1)/2. 51 ٢-7-٢7 
BROS((SIGC*SIGC*(X*QHAT-Z)* (X*QHAT-Z) ) *CCDEB1)/2. -SIGC/2. * 
*( X*QHAT-Z) *PBl 

EBO=(BR-BRQ) /QHAT 

RETURN 

END 
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I -Ὁ- Erud au XE مج و بو مه موجه دم‎ 
THIS IS A NONINTERACTIVE PROGRAM UTILIZING THE WATFIV COMPILER.THE * 
PROGRAM READS DATA WHICH WAS PREPARED BY THE INVENTORY PROGRAM e 
(APPENDIX A),AND THEN CALCULATES ISOCOST PARAMETERS FOR THE VARIOUS* 
CONBINATION OF C AND L WHICH HAVE THE SAME TOTAL COST. THE PROGRAM * 
IS ACTIVATED BY AN EXEC FILE . x 


ας κ em KKK KKK KKK KKK KK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK 


********** VARIABLE DEFINITIONS ***»******xxx*X 
S SLOPE OF THE ISOCOST LINE 

ANOM = NOMINATOR FOR CALCULATING A 

۳۱ INTERCEPTION POINT OF THE ISOCOST LINE 
BNOM - NOMINATOR FOR CALCULATING B 

SS= UNIT PRICE 

ee = C*C 

CL = C*L 

DENO = DENOMINATOR FOR CALCULATING A AND B 
J = LOOP COUNTER 

K = NUMBER OF TERMS IN THE REGRESSION 

KK = CURRENT TOTAL COST 

L = LEAD TIME 


BG = L*L 

LK = NUMBER OF C UNITS IN TC RANGE 

R = CORELATION COOFICIENT 

RDEN = DENOMINATOR FOR CALCULATING R 

RNOM = NOMINATOR FOR CALCULATING R 

SUMA = ACCUMULATIVE SUM OF A VALUES 

SUMB = ACCUMULATIVE SUM OF B VALUES 

SUMC = ACCUMULATIVE SUM OF C VALUES 

seme — ACCUMULATIVE SUM OF THE PRODUCT C*C 
NCE = ACCUMULATIVE SUM OF THE PRODUCT C*L 


SUML = ACCUMULATIVE SUM OF L VALUES 
ο ΕΕ = ACCUMULATIVE SUM OF THE PRODUCT L*L 
ne TOTAL COST 


Ke kkk KKK Ke VARIABLE DECLARATION *********** 
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INTEGER SUMC,SUMCC 
INTEGER C(1999),TC(1999),CC(1999) 
REAL A,B,L(1999),5L(1999),CL(1999), SUML, SUMLL, SUMCL 
C *** READ INPUT DATA *** 
DO 40 J=1,960 
MEAD C1 C(I). L(JJ, TCC) 
41 FORMAT(T3,I4,T8,F5.2,T15,I5) 
40 CONTINUE 
KK=TC(1) 
۱۰۰-0 
K0 
20 SUML=O 
SUMC=O 
SUMIC=0 
SUMLL=0 


SUMCL=0 

K=0 

LOOP OVER ALL TC VALUES *** 
DO 10 J=1,960 

τη (τοι) NE. KK common ma 
LINEAR REGRESSION COMPUTATIONS *** 
SUML=SUML+L(J) 
SUMC=SUMC+C(J) 

ασ 7٦ 

WU CODI δα. 

CEC LT) 

CE=C( a) 

SUMCL=SUMCL+CL(J) 
SUMCC=SUMCC+CC(J) 
SUMLL=SUMLL+LL( J) 


K=K+1 

IE ας NE T CO TON 
ClsC(J) 
CONTINUE 
CONTINUE 


BUCSUGI OE) 22. 
IF(K LE 1 21 
DENO=K* SUMCC-SUMC* SUMC 
IF (DENO.EQ.O)GO TO 21 
BNOM=K * SUMCL-SUMC*SUML 
ANOM=SUMCC* SUML-SUMC* SUMCL 
BB=BNOM/DENO 
BBB=BB*100. 
B=AINT(BBB)/100. 
A=ANOM/DENO 
RNOM=K* SUMCL- SUMC* SUML 
RDEN- ( KKSUMCC-SUMC*SUMC ) * ( K*SUMLL-SUML*SUML) 
R-RNOM/SQRT(RDEN) 
PRINT 23,KK,A,B,K,R,AVC 
FORMAT (2X, I5,2X,F7.2,2X,E7.2,2X, 12, 2K E ee 
IE(AVC NE CAC) COTTO IQ) 
CALCULATIONS OF A AS FUNCTION OF C *** 
SUMA=SUMA+A 
SUMB=SUMB+B 
LK=LK+1 
οσο πο ον 
و ئ۲‎ ٤ ۶ 
AA=SUMA / LK 
AB=SUMB / LK 
PRINT,CAC,AA,AB 
CAC=AVC 
SUMA=A 
SUMB=B 
ΤΕ 
KK=KK+10 
IF(KK OCIO TOUJ20))60 T0950 
50, ΤΟ. 28 
CONTINUE 
STOP 
END 


C *** 
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